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contains  two  features  that  will  be  of  aid  to  the  programer.  These  are  sample  □ 
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programs,  written  in  BASIC  for  a  Casio  programmable  calculator,  and  numerical 
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Cartographers  involved  in  national  mapping  t.i  rams  frequently  need  to 
transform  geographic  coordinates,  by  way  of  a  specific  projection  system,  to 
their  cooresponding  rectangular  coordinates.  They  may  want  to  plot  control  points 
or  construct  graticules  of  meridians  and  of  latitude. Until  recently  this  was  a 
time  consuming  chore.  Now,  with  the  proliferation  of  hand  held  microcomputers 
and  programmable  calculators,  this  once  tedious  job  is  performed  quickly  and 
accurately .  In  the  time  it  took  to  transform  one  set  of  geographic  coordinates 
into  rectangular  coordinates  using  projection  tables  we  now  can  perform  many. 

But  there  is  one  problem  that  is  holding  back  the  greater  use  of  hand  held 
computers.  T'’is  problem  is  the  difficulty  in  locating  concise,  easy  to  program 
transformation  formulas. 

The  purpose  of  this  article  is  to  provide  Cartographers  with  the  transform¬ 
ation  equations  for  the  most  widely  used  grid  system,  the  Universal  Transverse 
Mercator  grid  system  (UTM  grid). 

The  equations  perform  forward  transformations;  geographic  coordinates  to 
UTM  grid  coordinates,  and  inverse  transformations;  UTM  grid  coordinates  to 
geographic  coordinates.  These  equations  are  notable  for  the  ease  of  which  they 
car.  be  programed. 

The  article  begins  with  a  short  history  of  the  development  of  the  UTM  grid 
system  and  then  states  its  basic  characteristics.  Following  this  the  forward 
and  inverse  transformations  are  presented.  The  last  section  of  this  article  ___ 
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The  history  of  the  development  of  the  *JTf<*  grid  system 
The  Transverse  Mercator  projection  is  used  in  the  UTM  grid  system.  The  Transverse 
Mercator  projection  evolved  from  the  regular  Mercator  projection,  named  after 
the  Flemish  cartographer  Gerhardus  Mercator  (Gerhard  Kremer,  1512  -  1594). 

The  development  of  the  Mercator  projection  can  be  represented  geometrically  by 
viewing  a  cylinder  that  is  wrapped  around  the  globe  and  tangent  to  the  Equator. 

(see  figure  1).  A  light  source  approximately  3/4  of  the  way  back  along  the 
diameter  can  be  seen  as  projecting  tne  earths  features  out  on  to  the  cylinder. 

The  cylinder  is  then  cut  and  laid  out  flat  to  form  the  map.  The  feature  that 
is  interesting  to  us  here  is  that  there  is  very  little  scale  distortion  close 
to  the  standard  line,  the  Equator. 

The  Transverse  Mercator  was  developed  to  take  advantage  of  this  feature.  If 
you  take  the  cylinder  and  rotate  it  90“  you  effectively  trace  a  meridian  the 
standard  line  (see  figure  1).  This  makes  the  Transverse  Mercator  useful  in 
mapping  narrow,  North-South  orientated  features. 

The  spherical  form  of  the  Transverse  Mercator  projection  was  first  developed 
by  Johann  Heinrich  Lambert  (1728  -  1777),  an  Alsatian  mathematician,  astronomer 
physicist  and  philosopher.  His  work,  Beitr’ege,  published  in  1772  contained  six 
new  map  projections  and  included  the  Transverse  Mercator  and  the  Lambert  Conformal 
Conic.  It  is  interesting  to  note  that  the  Lambert  Conformal  Conic  and  a 
modified  -ell ipsoical  form  of  the  Transverse  Mercator,  tne  UTM,  are  the  two 
most  widely  used  projections  for  large  and  medium  scale  mapping  in  Latin  America 
today. 

After  Lambert,  the  Transverse  Mercator  was  worked  on  by  the  great  German 
mathematician,  Johann  Karl  Friedrich  Gauss  (1777  -  1855).  It  wasn't  until  1912  and 
1919  that  L.  Kruger  published  his  studies  that  for  the  first  time  stated  the 
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Transverse  Mercator  ellipsoidal  formulas.  Because  of  the  work  of  Gauss  and  Kruger 
the  Transverse  Mercator  is  sometimes  referredto  as  the  Gauss-Kruger  projection. 

Tne  L’TM  grid  system,  which  is  based  on  the  ellipsiodal  form  of  the  Transverse 
Mercator  projection,  grew  out  of  the  military  need  for  a  worldwide  plane 
coordinate  system.  The  ideal  universal  plane  coordinate  system  was  seen  as 
containing  the  following  properties;  conformality,  continuity  over  largo  areas, 
a  minimum  number  of  zones,  minimal  scale  errors,  a  unique  referencing  svsLom, 
uniform,  transformation  formulas  and  a  maximum  meridional  convergence  of  five 
degrees.  Based  or.  these  criteria,  in  the  late  forties  the  U.S.  Army  Map  Service 
devised  the  UTM  grid.  Today  the  UTM  grid  system  and  its  companion  system, the 
Universal  Polar  Stereographic  (UPS)  ,  for  use  in  the  polar  regions, are  used 
extensively  in  the  Americas,  Europe  and  Russia  for  surveying,  mapping  and 
military  purposes. 

Basic  characteristics  of  the  UTM  grid  system 

The  UTM  grid  system  is  made  up  of  sixty  zones  each  covering  six  degrees  >  f 
longitude.  The  zones  are  numbered  one  through  sixty  starting  at  180°  west  from 
Creenwici.  and  increasing  to  the  east.  The  zones  extend  in  latitude  from  t.i°N 
tc  SC'S  (sec  figure  2). 

Each  UTM  grid  zone  is  bisected  by  a  central  meridian  which  becomes  the  V 
axis  origin  of  the  projection.  The  central  meridian  is  assigned  a  value  of 
50'  ,uGC  meters  in  order  to  avoid  the  use  of  negative  numbers.  Along  the  central 
meridian  a  scale  factor  of  0.9996  is  employed.  Because  of  this  the  central 
meridian  is  not  a  standard  line.  The  effect  of  using  a  scale  factor  of  less  then 
one  is  to  create  two  standard  lines,  in  this  case  standard  meridians,  spaced 
1  SC  km.  east  and  west  of  the  central  meridian.  By  having  two  standard  meridians 
the  scale  retention  characteristics  of  the  UTM  are  greatly  improved. 
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(see  figures  3  and  4) . 

The  X  axis  origin  is  the  Equator.  When  working  in  the  northern 
hemisphere  the  Equator  is  assigned  a  value  of  zero.  In  the 
southern  hemisphere  the  Equator  is  given  a  false  value  of 
10,000,000  meters.  The  false  value  is  appli..-  to  avoid  the 
use  of  negative  numbers. 

The  UTK  Eouations 

The  following  formulas  are  taken  from  the  publication, 

"Map  Projections  used  by  the  U.S.  Geological  Survey"  second 
edition,  1984.  They  are  presented  here  in  a  modified  form. 

Symbol s 

a=  Semimajor  axis  of  the  reference  ellipsoid 
b=  Semiminor  axis  of  the  reference  ellipsoid 
e=  Eccentricity  =  (l-fc2/a2)j. 

x=  Rectangular  coordinate:  the  distance  to  the  right  of  the 
vertical  line  (the  y  axis)  passing  through  the  origin  of  the 
projection,  the  central  meridian. 

y=  Rectangular  coordinate:  the  distance  above  the  horizontal 
line  (the  x  axis)  passing  through  the  origin  of  the  projection, 
the  Equator. 

k=  The  relative  scale  factor  along  a  parallel  of  latitude, 
ke  =  Tne  scale  factor  of  the  central  meridian.  K<.=  0.9996. 
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0  =  The  latitude  of  the  point.  For  a  latitude  south  of  the  Equator 
apply  a  minus  sign. 

0O=  The  latitude  of  the  origin.  0°  in  the  UTM  system. 

"X =  The  longitude  of  the  point.  For  longitude  west  of  Geenwich 
apply  a  minus  sign. 

=  The  longitude  of  the  origin,  the  central  meridian.  For  longitude 
west  of  Greenwich,  apply  a  minus  sign. 


Northing  =  y,  in  the  northern  hemisphere 

Northing  s  10,000,000  +  y,  in  the  southern  hemisphere 


x=k0N[A  +  (1  -  T  +  C)A3/6+{5  -  18T  +  T2  +  72C  -  58  e'2)A5/120] 

y=k0£f-HN  tan  0  [A2/2+(5  -  T  +  9C  +  4C2)  A4/24+(6l  -  58T  +  T2  +  600C-330e'2) 
A6/720]> 

k=kon+(l+C)A2/2+(5  -  4T  +  42C  +  13C2  -  28e'2)  A4/24+(61  -  148T  +  l6T2)A6/720 
Where  ko=0 . 9996 


e  2=e2/(l-e2) 

N=a/(l-e2sin2  0)* 

T=tan20 

C=e'2cos20 

A=cos  0  (  X  -  X  o) .  with  ^  and  Xo  in  radians. 
M=a[(l-e2/4-3e4/64-5e6/256-. . . )  0-(3e2/8+3e4/32+45e6/lO24+. . . ) 
sin20+(15e4/256+45e6/lO24+. . . )  sin40-(35e6/3O72+. . . )sin60+. . .] 


With  -0  in  radians . 


THE  INVERSE  EQUATIONS  FOR  THE  UTH  GRID 


0  =  0]-(N1  tan  01/R1)[D2/2-(5+3T1  +  10^-40^  -  9e'2)  D4/24  +  (61  +  90T-J  + 

298C i  +  45T-J2  -  252e  ‘"-3C]2)D^/720]  ,  wnere  0  is  in  radians.  To  convert  to  degrees 
mul  tiply  0  Dy  180/TT  . 

Ac  +  [D-  (1  +  2T -j  +  C-|)  D3 / 6  +  (5-2C1+28T1  -  30^  +  8e'2  +  247-, 2 )' 

D^/ 120 j / COS  0i  ,  where  A  is  in  radians.  To  convert  to  degrees  multiply  A  by  1  OOAJT 
Where  0i  is  tne  foot  point  latitude  of  the  perpendicular  from  the 
central  meridian. 

01  =  m  +  (3e1/2-27e13/32+...)  sin  2u  +  ( 21  e-, 2/l 6  -  55e-,4/32-. . . ) 

sin  4u  +  ( 1 51  e-,  /96+...)  sin  6u  wnere  0  is  computed  in  radians  and  then 

converted  to  aegrees  by  multiplying  by  1 8O/77-"  . 

Where 

ei  =  [l-(l-e2)-]/[l+(l-e2)-] 

i x  =  M/[a(l-e2/4-3e4/64-5e^/256-. . . )]  ,  where  tx  is  in  radians. 

*  =  y/*0 

kc=  .9996 

With  0-  known,  tne  otner  terms  are  calculated  for  use  in  the  equations  0  and  X. 
e'-=e2/'l-e2) 

:•  =  e  ^ccs-01 
=  tar70« 


The  UTK  coordinate 
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transformation  procrams 

The  two  programs  presented  are  written  in  BASIC  for  the  Casio 
FX-801P  Programmable  Calculator.  They  perform  the  forward  .and 
inverse  transformations. 

The  programs  operate  interactively.  The  forward  transformation 
begins  by  inquiring  of  the  user  which  reference  ellipsoid  he  wishes 
to  use.  This  is  accomplished  by  the  interrogative  "A=”  and  "E 
C'JADRADO="  .  A  is  the  semimajor  axis  of  the  reference  elliosoic  and 
E  C'JADRADO  is  the  eccentricity  squared.  These  values  for  four 

reference  ellipsoids  common  to  Latin  America  are  given  in  table  1 . 

2 

A  and  e  need  to  be  entered  only  once.  The  program  will  skip  these 
inquiries  or.  the  second  point  transformation.  After  entering  the 
ellipsoid  constants  the  program  will  ask  for  the  "LAT  DEL  PUNTC"  and 
the  " LON  DEL  PUNTO the  latitude  and  longitude  are  entered  in  the 
following  manner.  To  enter  35°  22'  14.22”  K  you  press  25  then  execute, 
22  then  execute  and  14.22  and  execute.  To  enter  a  point  located 
in  tne  southern  hemisphere  or  west  of  Greenwich,  you  must  apply 
a  minus  sign.  75°  34'  2C‘  would  be  entered  as  -75  execute,  -34 
execute,  -20  execute.  Once  the  latitude  and  longitude  are  entered 
the  program  will  ask  for  the  "MEP.IDIANO  CENTRAL"  .  This  is  entered 


as  two  digi 

ft 

tn 

Per  example  £9°K  will  be  entered 

as 

-69 

execu 

te  . 

After  this 

the 

transformation  calculations  will 

be 

oer 

f  ernee 

ciTiC 

the  answer  displayed. 

The  display  sequence  is  as  follows;  *’K=”  and  the  Northing  will 


After  the  Eastmc  is  displayed  the  prograrr  ’./ill  begin  aaam  bv 
asking  for  the  latitude,  longitude  and  the  central  meridian  of  the 
next  point. 

The  inverse  transformation  program,  like  the  forward  program; 
starts  by  asking  for  ,,A="  and  "E  CUADRADO="  .  The  program,  ' then  asks 
for  the  Northing  as  "NORTE  CUTM=" ,  and  the  Easting  as  "ESTE  CUTK " . 
After  this  the  program,  will  ask  for  the  "MERIDIANO  CENTRA.!".  This 
is  entered  as  two  digits,  like  in  the  forward  program.  Finally, 
the  program,  will  ask  for  the  hemisphere  in  which  the  point  lies. 

The  statement  "HLMISFLKiU  N  U  b"  will  appear.  Enter  if  you  are 

m  the  northern  hemisphere  or  S  for  the  south.  After  the  trans¬ 
formations  are  complete'  the  results  will  be  displayed  in  the 
following  manner.  "LATITUD  ="  will  appear.  Press  the  continue  key 
tc  display  the  results.  Press  continue  again  and  "LONGITUD  = "  will 
appear.  Press  continue  to  see  the  results.  If  you  wish  to  trans¬ 
form.  more  points  press  continue  and  the  program  will  ask  you  for 
the  Northing,  Easting,  central  meridian  and  hemisphere  of  the  next 
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The  UTK  Forward  Proaram 


The  UTM  Inverse  Proaram 


I 
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Programs  written  by  Liam  P.  O'Brien 
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Nii"ior  i  r.i !  Exnmo]  e*. 


”TM  Forward  Equations 

Given:  Clarke  1866  ellipsoid:  a=6378206. 4m 


Central  meridian: 
I’TM  Zone  18: 

Central  meridian 
Scale  factor: 

Point : 


e"=0. 00676866 


?\0  =  15c  West  Ion'. 


k,  =0 . 9996 


P=40D30'  North  latitude 
?\=7  3°  30 '  West  longitude 


Find:  x,v,k 


e  ’  '=0. 00o7ti8nb/  (  1-0. 00b768hb)=0. 0068  I  48 
N  =6378306. 4/ ( 1 -0. 00676866  sin2  40 . 5° ) !;=b387 330 . 3m 
T  =  tan 2 40.  r>c=0. 7. ’64 6  38 
C  =C.00o8U8cos240.5c=0. 0039404 
A  =(cos40°30')*:[  (-'3. 5°)- (-75°)  In'/!  80c=0.  01  99074 

y.  *637  6206. 4  (  ( 1-0. 0067  6866/4- 3x0. 006768662/64-5:-:0. 0067  6866-/ 2  36  ):■(• -0. 2C 
( 3x0 . 00676866 / 8+3x0 . 0067 68662 / 32+45x0 . 0067 6866** / 1 024 )  y.s  i  n  ( 2  >:4C .  5  0  >- 
C5:sC.006768662/256-ri5x0.006768664/1024)vSin(4x4C.3o  .-(23>:C.006‘,666t"/30'2): 
sin (6x40. 5°) Iff/ 180*  *-,484,837. 67m 

x  =0. 999 6x63873 30. 5x[0. 019907 4+ (1-0. 7294538+0. 00 3 9404) xC . 0; 9907- ' /t— 

(5-1 SxC. 7294538+0. 72945382+72x0.  0039404-58x0. 00681 48) xC .  C199C743.- ;  20 ; 

=1 27 , 1 06 . 5m 

v  =0. 9996x{-484£37 . 7-0+6387  330. 5x0. £5408G7x [ 0.01 99074  "72+ (5-G . "294 ' 3r+9x 
0. 0039404+4x0. 0C 39404 2 )x0. 0! 99074^ /24+( 6 1-58x0. 72945538+C .724-338 
600x0.0039404-330x0. 00681 48 )x0. 0199074^/720]} 

=4 , 484 ,124.4m 

k  =0.  999ox[  1  +  O+0.  00681 48cos240. 5°  )xl  27 1 06.  5'/  (2x0 . 99962x6  387  330 . 5  2 )  ) 

=0.9997989 

Northing  =4 , 484 , 1 24 . 4m 

Easting  =500 , 000+1 27 , 1 06. 5=627  ,! 06.  5m 


...  >  -  >  ji  r-'.r.j*  r n  ‘j-  -- *•  *  . *  X"  ->  .  -  ■  w  "j-  - 


UTM  Inverse  Equations 
Given:  Clarke  1866  ellipsoid: 


CM 


Central  meridian: 
UTM  Zone  18: 

Central  meridian 
Scale  factor: 

Point : 


a=6378206. 4m 
e2=0. 00676866 

Ae=750  West  longitude 

ko=0.9996 

Northing  *  4484124.4  m 
Easting  *  627106.5  m 


Find: 

>;  =  627106.5  -  500000  =  127106.5  m 
y  *  448m  124 . 4  m 

e  '  2  *0.00676866/(1  -0. 00676866 )=0. 0068148 
M  =0+4484124. 4/0. 9996=4485918. Sm 

e,  =[ •-( 1-0. 00676866)^] /[ l  +  ( 1-0. 0067  6866)^] =0.001 69791  6 

M  =4485 9 16. 8/ [ 6378206. 4x( 1-0. 0067 6866/4-3x0. 0067 68662/64-5x0.006768663/256) ] 
=0.7045135  radian 

01  *=C. 7045 135+(3x0. 001697916/2-27x0. 0016979163/32)sin(2x0. 7045135)+ 

(21x0. 001 69791 62/16-55x0. 001697916^* /32)sin (4x0. 70451 35)+( 151x0. 001 6979164/96) 
sin (6x0. 70451 35 )=C. 7070283  radian=C. 7070283x180''  ,*=40 . 5097362° 

C,  =C. 0068 14&cos-40. 5097362 '=0.0039393 
7-  =tar24C.5097362c=C. 7299560 

K,  =6378206.^/(1-0. 00676866s ir. "AO. 5097362"  1^6387 334.2m 

R:  =6378206. 4x(l -0 . 00676866)  /  ( 1-C.  00676866sir.240. 5097362°)2/,2=6 , 362 ,271.4m 

D  =127 106. 5/(6387334. 2x0. 9996) =0.0199 07 7 

P  =40. 5097362 ‘-(6387334. 2x0. 8543746/636227 1.4)x[ 0.01 990772/2- (5+3x0.7299560+10 
xC. 0039395-4x0. 0039393:-9x0. 00661 48) xO. 01 990774/24+(61+90x0. 7299560+2 98 
xC. 0039393+45x0. 729?56C2-252x0. 00681 48- 3x0. 00393932)x0.01990776/720xl80V 
=40. 500000CC=40°3C’ North  latitude 


MM 


A  =-? 01 99077- (1+2x0. 7299560+0. 0039393)x0. 01 99077 3/6+(5-2x0. 0039393+ 
28x0. 7299560-3x0. 00 393 932+£x0. 0068148+24x0. 7 2995 602)x0. 01 990775/ 1 20] /cos 
40.5097362eJxl8CoA'  —  75e+1.5000000t=-72.5c=73o30’  West  longitude 


Table  1 


Constants  for  Reference  Elliosoids  Common  to  Latin  America 


Reference  Ellipsoid  Semi  major 

Clarke  1866  6378206.4 

Internacional  6378388 

South  American  1969  6378160 

i;GS  72  6378135 


Axis  a  Excentricity  Squared  e^ 
0.0067686580 
0.0067226700 
0.0066945419 
0.0066943178 


t.n 
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Figure  1 


Regular  Cylindrical 


Standard  line 


Transverse  Cylindrical 


The  cylindrical  developable  surface  as  seen  in  its  regular 
and  transverse  forms. 


BMC 


The  Central  Meridian. 

The  Standard  Lines  formed  by  the  intersections 
of  the  cylinder  and  the  spheroid. 


A  C  D 


r*. 


